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O ■ ABSTRACT 

<N ■ 

D ■ We constrain the X-ray emission properties of Lyman break galaxies (LBGs) at z ~ 3-6 

using the ~ 2 Ms Chandra Deep Field-North and fh 1 Ms Chandra Deep Field-South. Large 
samples of LBGs were discovered using HST as part of the Great Observatories Origins Deep 
Survey (GOODS). Deep optical and X-ray imaging over the GOODS fields have allowed us 
to place the most significant constraints on the X-ray properties of LBGs to date. Mean X- 
ray properties of 449, 1734, 629, and 247 LBGs with z ~ 3, 4, 5, and 6, respectively, were 
determined using stacking techniques. When stacked, we detect X-ray emission from LBGs at 
Q\ • z ~ 3 (~ la) and from an optically bright subset (brightest 25%) of LBGs at z ~ 4 (~ 3a); the 

latter is the highest redshift detection yet for "normal" galaxies in the X-ray band. The effective 
exposure times for these stacked observations are « 0.7 and 0.5 Gs, respectively. The derived 
"c^ average rest-frame 2.0-8.0 keV luminosities are 1.5 and 1.4 x 10 41 erg s" 1 , respectively. X- 

ray emission from these LBGs is likely due to high mass X-ray binaries (HMXBs) and Type II 
supernovae; the corresponding star formation rates are k, 85-240 M yr" 1 . The X-ray to fi-band 
' mean luminosity ratio (L X /L B ) at z, ~ 3 is somewhat elevated with respect to that measured for 

starburst galaxies in the local Universe (significance ~ 3a). When stacking full samples of 
LBGs at z ~ 4, 5, and 6 we do not obtain significant detections (< 3a) and derive rest-frame 
^ . 2.0-8.0 keV luminosity upper limits (3a) of 0.9, 2.8, and 7.1 x 10 41 erg s , respectively. 

These upper limits constrain any widespread AGN activity in these objects to be modest at best. 
Furthermore, we find that ~ 0.5% of our LBGs from z ~ 3-6 are detected individually in the X- 
ray band. These LBGs have spectral shapes and luminosities characteristic of moderate-power 
AGN (e.g., Seyfert galaxies and quasars). 

Subject headings: cosmology: observations — cosmology: surveys — X-rays: galaxies — 
X-rays: general 
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1. Introduction 

Determination of the basic X-ray properties of "normal" (i.e., not hosting a luminous AGN) galaxies 
out to high redshift has been motivated in part by Ghosh & White (2001, hereafter GW01) who developed 
a model of the evolution of X-ray luminosity with redshift. GW01 predict that, globally, X-ray emission 
from normal galaxies should peak at z ~ 1.5-3 due to a maximum in the global star formation rate (SFR) 
at z ~ 2.5-3.5 (e.g., Blain et al. 1999). Heightened X-ray emission at these redshifts is expected due to an 
increase in the global population of low mass X-ray binaries (LMXBs). LMXBs evolve on timescales of 
~ 1 Gyr, and therefore their X-ray signatures should lag behind more immediate tracers of star formation. 

Stacking analyses and deep X-ray surveys with the Chandra X-ray Observatory (hereafter Chandra) 
have enabled testing of these predictions out to cosmologically interesting distances. Stacking techniques 
generally involve the addition of X-ray counts at known positions of galaxies to yield average X-ray detec- 
tions where individual sources lie below the detection threshold. For example, Hornschemeier et al. (2002, 
hereafter H02) used stacking with a « 1 Ms exposure of the Chandra Deep Field-North (CDF-N) to inves- 
tigate the evolution of the X-ray luminosities of z = 0.4-1.5 spiral galaxies. This analysis shows suggestive 
evidence for evolution in X-ray luminosity with redshift at a level somewhat lower than that predicted by 
GW01. H02 also found that the X-ray to fi-band luminosity ratio, Lx/Lb, is elevated for z ~ 1 spirals with 
rest-frame Lb ~ L* B as compared with galaxies in the local Universe. Lx/Lb has been shown to be linked to 
star formation activity. This is supported by the existence of an empirical correlation between Lx/Lb and 
^60/im/^i00/im (Fabbiano & Shapley 2002). Here, the quantity LeO/um/^ioOjum is a measure of far infrared 
color temperature, which has been shown to increase with star formation activity (e.g., Helou, Ryter, & 
Soifer 1991). More recent stacking analyses of large samples of galaxies in the m 2 Ms CDF-N and « 1 Ms 
Chandra Deep Field-South (CDF-S) show that the X-ray properties of normal galaxies of all morphological 
types evolve similarly from z = 0.4-1.5 (Wu et al. 2004). 

At higher redshifts (z ~ 2-6) stacking analyses have been performed using individually undetected nor- 
mal galaxies (e.g., Brandt et al. 2001, hereafter B01; Nandra et al. 2002, hereafter N02; Malhotra et al. 2003; 
Bremer et al. 2004; Moustakas & Immler 2004 Reddy & Steidel 2004; Wang et. al. 2004). At z ~ 3 these 
galaxies often have mean X-ray luminosities and Lx/Lb comparable to bright starburst galaxies in the local 
Universe. The heightened X-ray emission from these galaxies is probably due to newly formed high mass 
X-ray binaries (HMXBs) and young supernova remnants. HMXBs have much shorter formation timescales 
than LMXBs and are therefore more immediate tracers of cosmic star formation history. 

Recently, large samples of galaxies have been identified at z ~ 3, 4, 5, and 6 as part of the Great 
Observatories Origins Deep Survey (GOODS). These galaxies were isolated via the Lyman break technique 
(e.g., Steidel et al. 1995; Madau et al. 1996; Steidel et al. 1999), which estimates the redshift of a galaxy 
based on its "dropout" bandpass (see § 2.1). Lyman break galaxies (LBGs) generally show rest-frame 
UV/optical characteristics similar to those of local starburst galaxies — weak or absent Lyman a emission, 
P Cygni features from C IV and other lines, and strong interstellar UV absorption lines such as Si II, O I, C II, 
Si IV, and Al II (see e.g., Steidel et al. 1996; Shapley et al. 2003). The GOODS covers « 316 arcmin 2 of sky 
in two fields, GOODS-North (GOODS-N) and GOODS-South (GOODS-S) (Giavalisco et al. 2004a). Both 
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fields are subregions of the Chandra Deep Fields for which Chandra has acquired an « 2 Ms observation 
of the CDF-N (Alexander et al. 2003a, hereafter A03) and an ps 1 Ms observation of the Chandra Deep 
Field-South (CDF-S; Giacconi et al. 2002; A03). 

In this paper, we use stacking techniques to constrain the X-ray properties of LBGs in the GOODS 
fields. We improve on results from B01 and N02 by increasing the number of galaxies at z ~ 3 from 24 
(B01) and 148 (N02) to 449, and we add to this samples of 1734, 629, and 247 galaxies at z ~ 4, 5, and 6, 
respectively. 

The Galactic column densities are 1.3 x 10 cm" " 2 (Lockman 2003) and 8.0 x 10 19 cm -2 (Stark et al. 
1992) for the CDF-N and CDF-S, respectively. H = 70 km s _1 Mpc -1 , Q m = 0.3, and S7 A = 0.7 are adopted 
throughout this paper (Spergel et al. 2003). Coordinates are J2000.0. 



The Lyman break technique has been utilized in the GOODS regions, where deep mosaic images 
have been taken using the HST Advanced Camera for Surveys (ACS) with bandpasses B435, V^oe, his, 
and Z850 (Giavalisco et al. 2004b, hereafter G04); the full-depth, five-epoch ACS images have been used 
here. Additional ground-based photometry has been obtained over the GOODS fields in the U band at 
KPNO (Capak et al. 2004, GOODS-N) and CTIO (GOODS-S). Identification of these LBGs is based on the 
following color equations: 



2. Analysis 



2.1. Samples 



U -dropouts: 



(t/-fl 4 50)>0.75 + 0.5(fl 4 50-Z850) A 



(U-B 450 )>0.9 A (fl 4 50-Z85o)<4.0 



(1) 



B435 -dropouts: 



(B450 - V 606 ) > 1 .2 + 1 .4 x (V 606 - Z850) A 



(#450-^606) > 1-2 A (V 606 -Z850)< 1-2 



(2) 



V606-dropouts: 



[(V606-«775) > 1-5 + 0.9 X (J775-Z850)] V 



KV606-l775)>2.0] A (V 606 -i775)>l-2 A 



(«775-Z850)<l-3 A S/N (fi 4 50)<2.0 



(3) 



1775 -dropouts: 



(1775 -Z850) > 1-3 A 



S/N(B 45 o)<2.0 A S/N(V 606 )<2.0 



(4) 
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Here the symbols A and V correspond to the logical operators AND and OR, respectively. Monte Carlo 
simulations of the color-selection procedure (e.g., G04) indicate that these color criteria select galaxies with 
mean redshifts and la redshift ranges of z = 3.01 ± 0.24, 3.78 ± 0.34, 4.92 ± 0.33, and 5.74 ± 0.36 for 
U-, B435-, V606-, and ?775 -dropouts, respectively. Dropout lists created in G04 were inspected and filtered to 
minimize the number of interlopers and contaminants in the samples. The interloper fraction is predicted to 
be w 10% for U -dropouts and Z^-dropouts. This fraction increases to become significant for V<so6 -dropouts 
and 1775 -dropouts (« 20% and 40%, respectively; see Dickinson et al. 2004; G04) 

The LBGs used in our stacking analyses have average redshifts and Z850 magnitudes as shown in Fig- 
ure 1 ; for comparison the Z850 magnitudes of several other stacking analyses in this redshift range are also 
plotted. Galaxies in this survey probe relatively faint optical fluxes over a wider range of redshifts (z ~ 3-6) 
than those used for previous X-ray stacking analyses. The corresponding lookback times for U-, B435-, 
V606-, and /775-dropouts are 11.5, 12.1, 12.5, and 12.7 Gyr, respectively (see Hogg 2000). 

The stacking procedure used in our analyses was similar to that of B01 and N02; the intent was to obtain 
average X-ray properties (e.g., luminosities and spectral shapes). At each LBG position, photon counts and 
effective exposure times (from exposure maps, which include the effects of vignetting) were extracted from 
circular apertures and summed to give a total number of counts and a total effective exposure time. We 
excluded LBGs located within 10" of individually detected X-ray sources in A03; this avoids contamination 
by unrelated sources. This exclusion process led to the rejection of 99, 267, 127, and 31 U-, B435-, V^-, 
and /775-dropouts, respectively (« 18%, 13%, 17%, and 11% of the respective LBG populations). A small 
number of LBGs (18) were found to be within the positional errors of Chandra sources; their characteristics 
are given in § 3.1 and listed in Table 1. 

The background was estimated through Monte Carlo analysis using a background map (see § 4.2 of 
A03). Each LBG position was shifted randomly within 25" of the original position, and background counts 
were obtained for each new position and summed just as they were for the LBGs themselves. This procedure 
was repeated 10,000 times to obtain an accurate estimate of the local background and its dispersion. To 
verify our results were not biased by gradients in the local background, we also tried shifting in regions with 
radii of 10", 15", 20", 30", and 35" and found no material differences in our results. 

Our stacking procedure maximized the signal-to-noise ratio (hereafter, S/N) 1 by varying the aperture 
cell size for extracting X-ray counts and the off-axis angle (angle between a particular source and the av- 
erage CDF-N or CDF-S aim points) within which sources were included in stacking. This procedure was 
performed using U -dropouts because of their strong X-ray signal. We avoided stacking sources that were 
outside a certain off-axis angle due to the increased size of the point-spread function (PSF) and confusion 
with background in this region. In the maximization process we chose a specific aperture cell size and 
stacked sources within various off-axis angles (hereafter, inclusion radii) until S/N was maximized. Sources 
within the resulting inclusion radius were then stacked using variable aperture cell sizes. S/N was maximized 



'Here, S/N = (S-B)/B 05 , where S is the source plus background counts and B is the average background counts obtained from 
the Monte Carlo procedure. This approximation is accurate for \S-B\ <B and B > 20, which applies to all cases in our analyses. 
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for a specific aperture cell size, and the process was repeated iteratively until convergence. Figures 2a and 2b 
show the result of this process. We found that S/N peaks for an inclusion radius ?a 9'.0 and an aperture size 
of « 1"5 (as expected for small off-axis angles); these values are used throughout our analyses. 

Mean count rates $ were calculated by subtracting our best estimate of the local background counts 
B from the background plus source counts S and dividing by the total effective exposure time T [i.e., 
<J> = (S—B)/T]. To convert count rate to flux we used the Galactic column densities discussed in § 1 and 
adopted a power-law photon index of V = 2.0, appropriate for starburst galaxies (e.g., Kim, Fabbiano, & 
Trinchieri 1992; Ptak et al. 1999). Sources were stacked in both the GOODS-N and GOODS-S fields, so we 
used a statistically weighted Galactic column density. X-ray luminosities were calculated following Schmidt 
& Green (1986): 

L x = 4^ /x( i +z) r-2 ergs -i (5) 
where dt is the luminosity distance (cm), and fx is the X-ray flux (erg cm" 2 s" 1 ). 

2.2. Subsets 

The large number of LBGs in our sample enabled investigation of the contribution of specific subsets 
to the signal. We chose to divide the original dropout lists based on galaxy morphology (U -dropouts only) 
and rest-frame fi-band luminosity. We used the CAS morphology system (Conselice 2003, hereafter C03) 
to divide the [/-dropout list according to the observed-frame zg5()-band (rest-frame AA 2000-2500) concen- 
tration and asymmetry parameters. We note that the CAS parameters have been derived using rest-frame 
UV light and may have somewhat different physical interpretations than those quoted for the local Universe. 
Concentration (C) is proportional to the logarithm of the ratio of the radii containing 80% and 20% of the 
source flux: 

C = 5 x log (r 80 % /r 2 o%) (6) 

Galaxies with high C values (e.g., bulge dominated systems locally) are generally brightest in their central 
regions (Conselice, Gallagher, & Wyse 2002). Asymmetry (A) is defined for a galaxy by taking the absolute 
value of the difference in fluxes of the galaxy's image and its 180° rotated analog and dividing by the original 
image flux (see Conselice, Bershady, & Jangren 2000; C03, for details). Generally disk-dominated systems 
and galaxy mergers are observed to have high A values in the local Universe. Parameter values of A = 0. 1 
and C = 2.7 were chosen to divide the [/-dropout list roughly in half, resulting in lists of LBGs with A < 0.1, 
A>0.1,C<2.7, and C> 2.7. 

We also divided our original dropout lists into subsets based on rest-frame fi-band luminosity. Rest- 
frame B-band luminosities were calculated using a spectral energy distribution (SED) and applying K- 
corrections to the fluxes (e.g., Schneider, Gunn, & Hoessel 1983; Hogg 2002). ^-corrected ACS magnitudes 
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were used in these computations, thus the derived rest-frame B-band luminosities are strongly dependent on 
the SED choice. An SED for LBGs was used in calculating ^-corrections. This SED was derived from 
the Bruzual & Chariot (2003) solar metallicity models, with a Salpeter initial mass function, and aged by 
144 Myr at constant SFR. Optically bright and faint sublists were generated by dividing each dropout list 
(i.e., U-, B435-, V606-> and 1775 -dropouts) at its mean rest-frame B-band luminosity. The corresponding mean 
rest-frame fi-band luminosities used in these divisions were Lg s=s 5.3, 5.5, 4.1, and 5.5 x 10 43 erg s -1 
(Mb ~ -20.6, -20.5, -20.1, -20.6) for U-, B435-, V<so6-, and 2775 -dropouts, respectively. Because rest-frame 
fi-band luminosity distributions are asymmetric with skew tails toward high luminosity, splitting at the mean 
leads to two lists with unequal sizes. 



3. Results 
3.1. Detected Sources 

Over the GOODS fields seven [/-dropouts, eight #435 -dropouts, two V606 -dropouts, and one 1775- 
dropout were found to be coincident with Chandra sources within the positional errors (see Table l). 2 
Positional errors, photon counts, effective photon indices, and fluxes for these sources were determined in 
A03. The 2.0-8.0 keV X-ray luminosities were computed using eqn. 5. 3 We investigated the probability 
of obtaining a false match by shifting all of our LBG positions and checking to see if the new positions 
were coincident with Chandra sources. LBG positions were shifted by 10" in 16 different directions, and 
an average of « 2 detections were obtained for each direction. 

Of the 18 detected sources, we find that 14 have corresponding spectroscopic (Barger et al. 2003a; Cris- 
tiani et al. 2004; Szokoly et al. 2004) and/or photometric (Barger et al. 2003a; Mobasher et al. 2004; Zheng et 
al. 2004) redshifts. Of these 14 sources, four (J123621.0+621412, J123642.2+620612, J123701.6+621146, 
J033243.2-274914) have redshifts inconsistent with those determined by the Lyman break technique. We 
note that the redshifts derived for three of the four sources (J123621.0+621412, J123642.2+620612, and 
J123701. 6+621 146) are photometric redshifts, implying there are still uncertainties in their values. These 
may perhaps be dismissed as being "inconsistent." The spectroscopic redshift of z = 1.92 for [/-dropout 
J033243.2-274914 would therefore be the only remaining discrepant case. This LBG may be an interloper 
(see § 2.1 for probabilities) or falsely matched. Furthermore, the colors of some AGN may satisfy the 
color criteria used to select normal galaxies of differing redshifts. This may cause a discrepancy between 
spectroscopic and color-selected redshifts derived for these AGN. 

The detected sources are mostly AGN with moderate luminosities (i.e., Lx ~ 10 43 -10 445 erg s" 1 ), 



2 Note that the Lyman break technique is a statistical process for isolating galaxies at a given redshift. Due to the statistical 
nature of this process, we do not expect to recover all objects within the ACS flux limits that may reside at the redshifts under 
investigation here. 

3 For reference, the on-axis sensitivity limit for the 2 Ms CDF-N is ~ 2.5 x 1(T 17 erg crrT 2 s -1 in the soft band corresponding to 
rest-frame 2.0-8.0 keV luminosities of 2.0, 4.0, 6.8, and 10.4 x 10 42 erg s' 1 at z = 3, 4, 5, and 6, respectively (T = 2). 
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characteristic of luminous Seyfert type galaxies to quasars. Several of these sources have been character- 
ized in previous investigations. J123648. 0+620941 is the highest redshift quasar known (z = 5.186) in the 
CDF-N and -S fields (e.g., Barger et al. 2003b; Vignali et al. 2002). J033201.6-274327, J033218.2-275241, 
J033243.2-274914, J033244. 3-275251 (Szokoly et al. 2004), and J033209.4-274807 (Schreier et al. 2001) 
have been spectroscopically classified as broad-line AGN. J123701. 6+621 146 and J123647.9+620941 are 
both Very Red Objects (I-K > 4, VROs; Alexander et al. 2002a). J033229. 8-275 105 has been identi- 
fied as a putative type 2 QSO (Norman et al. 2002). Finally, three additional sources (J033239.7-274851, 
J033242.8-274702, and J033250.2-275252) have been cataloged as high-redshift QSO candidates (Cris- 
tiani et al. 2004). 

Following Alexander et al. (2002b) we assume detected sources with 2.0-8.0 keV Lx > 1.5 x 10 42 
erg s -1 or T < 1.0 are AGN. Only one source (J123642.2+620612) falls outside these criteria, leaving 17 
clear AGN. The derived AGN fractions for our [/-, B435-, Veo6-, and /775-dropouts are s=s 1.2%, 0.4%, 0.3%, 
and 0.4%, respectively; these fractions are lower than others found in previous investigations. For compar- 
ison, the AGN fraction derived from N02 is « 2.7%. In comparison to our B435-, Va)6-, and /775-dropouts, 
our [/-dropouts and the LBGs used by N02 are relatively luminous in the optical. The corresponding AGN 
fractions therefore suggest X-ray luminous AGN are generally optically luminous as well. 



3.2. Stacked Sources 

The results of our stacking analyses of individually undetected sources are found in Table 2. Hereafter, 
unless stated otherwise, we discuss the results in reference to the soft-band stacking analyses where we 
obtain significant detections; we do not obtain significant detections in the hard- or full-bands (this result 
is understandable due to the significantly lower background in the soft-band). The stacking procedure was 
repeated for the GOODS-N and GOODS-S fields individually; results from the two fields were statistically 
consistent. Generally, we detect the stacked emission from [/-dropouts (~ la) and all subsets generated 
therefrom. We do not obtain significant detections (S/N > 3a) for general samples of B435-, V(^(,-, and 
/775-dropouts, but we note a suggestive positive fluctuation (~ 2a) for our 5435-dropouts. We do, however, 
obtain a significant detection for the bright subset of fi435-dropouts (~ 3a). Figure 3a shows the distributions 
of counts obtained for the individual galaxies being stacked for both U -dropouts and the bright subset of 
B435 -dropouts. We obtained an average of 3.1 and 2.3 counts per cell for U- and bright B435 -dropouts, 
respectively. The summed numbers of source plus background counts were 1351 counts for [/-dropouts and 
926 for bright B435 -dropouts. In the 10,000 Monte Carlo trials where we shifted the aperture cells to random 
positions and summed background counts, we found that no trials produced > 1351 counts for [/-dropouts 
and only 6 trials produced > 926 counts for S435 -dropouts. Gaussian statistics predict « and 5 Monte 
Carlo trials should exceed the total source plus background counts for our analysis of U -dropouts and bright 
#435 -dropouts, respectively. Figure 3b shows the Monte Carlo distributions obtained when stacking random 
positions over 10,000 trials. In both cases the detection confidence levels are greater than 99.9%. Stacked 
and smoothed images are displayed in Figure 4. The images have effective exposure times of ss 0.7 and 0.5 
Gs (22 and 16 yr) for U- and bright i?435-dropouts, respectively. 
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Stacked Va)6- and /775-dropouts and subsets thereof produced no significant detections (i.e, S/N < 3d). 
We also attempted to merge the bright subsets of V6O6 -dropouts and /775-dropouts and failed to obtain detec- 
tions. X-ray emission constraints for these LBGs are tabulated in Table 2. 

4. Discussion 

Using the pa 2 Ms CDF-N plus pa 1 Ms CDF-S we have placed constraints on the X-ray properties 
of LBGs identified in the 316 arcmin 2 GOODS-N and -S fields. We have used X-ray stacking techniques 
on samples of 449, 1734, 629, and 247 individually undetected LBGs at z ~ 3, 4, 5, and 6, respectively 
(U-, B435-, V(,o6-, and /775-dropouts, respectively) to obtain their average X-ray properties. X-ray emission 
from LBGs is expected to be largely due to activity associated with star-forming processes (e.g., HMXBs 
and young supernova remnants) with a low-level contribution from low-luminosity AGN (LLAGN; e.g., 
Ho et al. 2001). With the intent to investigate normal/starburst galaxies we have rejected LBGs coinci- 
dent with known Chandra sources from our stacking analyses. Considering Chandra's superb sensitivity 
for detecting luminous AGN, we are confident that this rejection process has effectively removed strong 
accretion-dominated X-ray sources that would heavily contaminate our stacked signal. We have identified 
18 LBGs coincident with Chandra sources (see § 3.1). With the possible exception of one source we find 
these objects are moderately luminous AGN comprising pa 0.5% of the LBG population from z ~ 3-6. 

Using our stacking procedure (see § 2) we detect X-ray emission from LBGs at z ~ 3 and an optically 
bright subset of LBGs at z ~ 4. X-ray count rates derived from our soft-band detections and hard-band 
upper limits (3a) were used to constrain the spectral slopes of these LBGs. The 2.0-8.0 keV to 0.5-2.0 keV 
band ratio ($2.0-8.0 kev/$0.5-2.o kev) is calculated to be < 0.9 for [/-dropouts corresponding to an effective 
photon index lower limit of V > 0.8; a consistent, less tightly constrained photon index is also derived for 
our bright S435 -dropouts. This lower limit is consistent with our assumed Y = 2 photon index, which is 
expected for galaxies with high star formation activity. The derived average 2.0-8.0 keV X-ray luminosities 
for these LBGs are pa 1.5 and 1.4 x 10 41 erg s -1 for [/-dropouts and bright B435 -dropouts, respectively, a 
factor of pa 5-10 times higher than for typical starburst galaxies in the local Universe (e.g., M 82; Griffiths 
et al. 2000). Furthermore, our LBGs are a factor of pa 2 less X-ray luminous than those studied by B01 and 
N02. 

Assuming the X-ray luminosity functions of our individually undetected LBGs at high redshifts have 
similar functional forms to those of normal/starburst galaxies in the local Universe, we can estimate the 
expected median X-ray luminosities for these LBGs. This is important because, for a luminosity function 
with significant skewness, the median will differ from the mean. We have investigated this using the David, 
Jones, & Forman (1992, hereafter DJF92) sample of 71 normal/starburst galaxies in the local Universe and 
find the mean to median X-ray luminosity ratio £ x mean /£ x median ~ 5 The mean 2.0-8.0 keV X-ray lumi- 
nosity for these 71 galaxies is pa 6 x 10 40 erg s -1 , a factor of pa 2.5 lower than the mean X-ray luminosities 
of our [/-dropouts and optically bright subset of #435 -dropouts. If we assume the £ x mean /X x median ra tio for 
our LBGs is similar to that of the DJF92 sample, the corresponding median 2.0-8.0 keV X-ray luminosity 
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would be ss 3 x 10 40 erg s l . 

The strong U -dropout signal allows investigation of the contribution of specific subsets to the photon 
statistics (see § 2.2 for details and Table 2 for results). We find that the optically luminous U -dropouts con- 
tribute most of the X-ray flux (S/Nb r i g ht ~ 1.5 x S/Ndim). It is possible that the emission from rest-frame 
UV may be heavily obscured by dust as would be expected if LBGs were a "scaled-up" population of ul- 
traluminous infrared galaxies (ULIRGs; Goldader et al. 2002). The X-ray to B-band mean luminosity ratios 
for these LBGs suggests that this is likely not the case, and the intrinsic dust attenuation from LBGs is sim- 
ilar to that expected for star-forming galaxies (see Seibert et al. 2002). We find that relatively asymmetric 
sources (A > 0.1) dominate the photon counts. However, the rest-frame fi-band luminosities of these LBGs 
are somewhat elevated as compared with sources of low-asymmetry index (i.e., A < 0.1), and we there- 
fore offer no further interpretation. Furthermore, stacking of the two LBG subsets split by concentration 
index (C > 2.7 and C < 2.7) show that these two morphological subsets have similar X-ray properties. No 
additional information was extracted from our divisions on asymmetry and concentration. 

Recent investigations have found that hard X-ray emission largely associated with HMXBs can be 
used as a direct indicator of SFR (e.g., Bauer et al. 2002; Ranalli, Comastri, & Setti 2003; Grimm, Gilfanov, 
& Sunyaev 2003; Persic et al. 2004). If we assume the majority of the X-ray emission from our stacked 
detections of [/-dropouts and bright Z^-dropouts is unobscured X-ray emisssion from HMXBs, we expect 
the corresponding mean SFRs to be ss 85-240 Mq yr" 1 . These SFRs were derived using the linear SFR- 
L™™ eV relation given as equation (2) of Persic et al. (2004) assuming a correlation error of 20%. For 
comparison, the mean SFRs derived from the rest-frame 1500 A fluxes are estimated to be ss 65 and 10 M Q 
yr" 1 for our [/-dropouts and bright B435 -dropouts, respectively, without UV extinction. When a correction 
is made for dust extinction, the corresponding SFRs are ss 400 and 60 M Q yr" 1 (see G04). The reasonable 
agreement between the X-ray and UV extinction-corrected derived SFRs broadly supports our assumption 
that the X-ray emission is dominated by star-forming processes with little contribution needed from LLAGN. 

The X-ray luminosities for star-forming galaxies (such as LBGs) are potential tracers of global star 
formation history (Lilly et al. 1996; Madau et al. 1996; Blain et al. 1999). A global estimate of the 
SFR can be implicitly drawn from these LBGs by considering the mean X-ray luminosity per mean fi-band 
luminosity, Lx/Lb (mean quantities). 4 This ratio should be interpreted with some caution due to the possible 
non-linear relationship between Lx and Lg (i.e. Lx oc L#). The power-law index (a) for this relation has 
been reported to range from ss 1 (DJF92; Fabbiano & Trinchieri 1985) to as high as ss 1.5 (Shapley et al. 
2001) for galaxies in the local Universe. We find that Lx/Lb shows suggestive evidence for evolution with 
redshift as illustrated in Figure 5a. Here our [/-dropouts and bright fi435-dropouts are plotted and compared 
with the results for normal galaxies at z ~ (DJF92) and early-type spirals from z ~ 0.05-1.5 (Wu et al. 
2004). From these data we infer a peak in Lx/Lb at z ~ 1.5-3.0, consistent with predictions of GW01. 
These results, while based on larger samples of galaxies, are consistent with the near constancy of Lx/^uv 



4 In our stacking analyses we compute an average Lx for a given stack of source positions. We therefore refer to Lx/L B as 
meaning < Lx > / < L B >, where < Lx > and < L B > are mean quantities. When necessary for comparisons, data regarding 
Lx/Lb have been converted to this form. 
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at z ~ 1 and z ~ 3 reported by N02 for Balmer break galaxies and LBGs, respectively. In this description 
the global SFR is expected to peak at z ~ 2.5-3.5. During this epoch, X-ray emission is observed due 
to contributions from HMXBs and supernovae. Shortly after the UV luminous population evolves, X-ray 
emission will continue in the form of LMXBs resulting in an increase in Lx/Lb. 

When stacking LBGs at higher redshift, z ~ 5 and 6 (V606- and ins -dropouts, respectively), we do 
not obtain significant detections. The AGN contribution may be slightly more significant for these LBGs 
where the Chandra exposures are only capable of detecting individual objects with 2.0-8.0 keV luminosities 
> 10 43 erg s" 1 (implying many Seyfert type AGN might not be individually detected at these redshifts). In 
light of these limits we therefore place constraints on the average AGN content of these LBGs. The derived 
rest-frame 2.0-8.0 keV luminosity upper limits (3<r) for the V606- and /77s -dropouts are 2.8 and 7.1 x 10 41 
erg s . Such upper limits are characteristic of bright starburst or low-luminosity Seyfert type galaxies and 
are the tightest constraints yet to be placed on the X-ray properties of LBGs at z > 5. For comparison, X-ray 
analyses of 44 LBGs at z ~ 5 (Bremer et al. 2003) and 54 LBGs at z ~ 6 (Moustakas & Immler 2004) con- 
strain their average 2.0-8.0 keV X-ray luminosities to be less than 3 and 7 x 10 42 erg s -1 (3<r), respectively. 

Figure 5b shows the overall Lx/Lb vs Lx results for both individually detected and stacked LBGs. 
Stacking analyses show that the individually undetected LBGs at z ~ 3, 4, 5, and 6 have Lx/Lb ratios 
characteristic of local starburst galaxies. Our individually detected LBGs have much larger Lx/Lb ratios, 
which are expected for X-ray luminous AGN. 

We gratefully acknowledge support from STScI grant HST-GO-09425.26-A and NSF CAREER award 
AST-9983783 (BDL,WNB), the Royal Society (DMA), and NSF grant 03-07582 (DPS). We thank Stefan 
Immler and Wentao Wu for useful discussions and development of software. 
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Fig. 1. — Z850 magnitude vs. redshift for X-ray stacking analyses of LBGs. Filled (our LBGs) and open 
symbols (other investigations) represent the estimated mean zsso an d redshift with error bars showing the 
la spread for each quantity. The Brandt et al. (2001) measurements (+) represent the positions of individual 
LBGs with spectroscopic redshifts. 
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Fig. 2. — (a) The achieved signal-to-noise ratio (S/N) as a function of inclusion radius for stacking U- 
dropouts. Note that the number of stacked sources rises with increasing inclusion radius and is maximized 
at « 9'.0. (b) Signal-to-noise ratio (S/N) as a function of extraction radius. Here, the signal is strongest for 
an aperture radius of ss \'/5. This process of optimizing the S/N was achieved iteratively and appears here 
in convergence. 




Fig. 3. — (a) Histogram of the count distribution for [/-dropouts {upper panel) and fi435-dropouts {lower 
panel). The vertical dotted lines show the mean source plus background counts per aperture cell, and the 
vertical dashed lines show the mean background counts per aperture cell, (b) Results from Monte Carlo 
estimates of the background level for [/-dropouts {upper panel) and i?435-dropouts {lower panel). The 
detection level is represented with downward pointing arrows. 




Fig. 4. — Stacked soft-band images of 449 [/-dropouts (sa 0.7 Gs exposure) and 395 optically bright B435- 
dropouts 0.5 Gs exposure). Stacked emission from these LBGs is significantly detected in the soft-band 
with significances of ~ la ([/-dropouts) and ~ 3a (bright i?435-dropouts). The images are 15" x 15" (0"5 
pixel" 1 ) and were adaptively smoothed at 2.5a using the CIAO tool CSMOOTH. The faint "nebulosity" 
observed in the optically bright i?435-dropouts is attributed to smoothing over noise. The black circles are 
centered on our l"5 radius aperture cell that was used in the stacking analyses. 
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Fig. 5. — {a) The X-ray (0.5-4.5 keV) to fi-band luminosity ratio vs. redshift. Note the suggested maximum 
of Lx/Lb at z ~ 1.5-3 is consistent with predictions from the cosmic star formation history. Data used in 
this figure come from the following sources: z ~ (DJF92), z ~ 0.05-1.5 (Wu et al. 2004), z ~ 3 and 4 - 
our U -dropouts and B435 -dropouts, respectively, (b) X-ray to B-band luminosity ratio vs. X-ray luminosities 
of individually undetected stacked LBGs (+) and individually detected LBGs (□); solid squares represent 
broad-line AGN. The individually detected sources are AGN (with one possible exception) and have much 
larger Lx/Lb ratios than our stacked LBGs. X-ray luminosities were calculated assuming a power law with 
photon index T = 2, and fi-band luminosities were computed using ^-corrections with an LBG SED (see 
§ 2.2). 
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Note. — Col.(l): Chandra source name. Col. (2): Observed ACS dropout bandpass. Col. (3): Positional error as reported in A03. Col.(4): Angular offset between positions 
determined by Chandra and ACS. Col.(5): Spectroscopic ("s" superscript) or photometric ("p" superscript) redshift. Col.(6): Soft-band X-ray counts. Col.(7): Soft-band flux. 
Col.(8): Hard-band flux. Col.(9): Effective photon index. Here a value of 1.4 was assumed when photon statistics were too low to determine accurate values. Col.(10): Hard-band 
X-ray luminosity. Col.(ll): z 850 -band magnitude. 



Table 2. Stacking Results For Filtered Dropout Lists 
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Note. — The data apply to stacking analyses in the soft band. Col.(l): Description of the LBG sample stacked. "Bright" and "Dim" subsets were created by splitting the "General" lists at 
rest-frame B-band luminosities f» 5.3 and 5.5 X 10 43 erg s _1 (Mb ~ -20.6 and -20.5) for [/-dropouts and R^-dropouts, respectively. Col.(2): Number of sources being stacked. Col.(3): X-ray 
source counts obtained from stacking. Col.(4): Mean background counts obtained from Monte Carlo simulations. Col.(5): Poisson error of the background B° 5 . Col.(6): Signal-to-noise ratio 
(S-B)/B 0,5 . Col.(7): Total, stacked effective exposure time. Col. (8): Number of Monte Carlo trials (out of 10,000) that produced a background estimate > S. Col.(9): X-ray count rate. Col.(10): 
X-ray (0.5-2.0 keV) flux. Col.(ll): 2-8 keV rest-frame luminosity. Col. (12): Rest-frame B-band luminosity. Col.(12): Logarithm of X-ray to B-band luminosity ratio. 



